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Abstract 

We present the first analysis of liarmonic generation data where 
the full potential of the generalized nonlinear Kramers-Kronig (K- 
K) relations and sum rules is exploited. We consider two published 
sets of wide spectral range experimental data of third harmonic gen- 
eration susceptibility on different polymers, the polysilane (frequency 



1 



range: 0.4 — 2.4 eV) and the polythiophene (frequency range: 0.5 — 2.0 
eV). We show that, without extending the data outside their range 
with the assumption of an a-priori asymptotic behavior, indepen- 
dent truncated dispersion relations connect the real and imaginary 
part of the moments of the third harmonic generation susceptibility 
a;^"x''^^ (3wj ^1 ^)-, with a ranging from to 3, in agreement with the 
theory, while there is no convergence if we choose a = 4. We repeat 
the same analysis for uj'^°^[x^'^\'iuj;uj,uj,ujy\^ and show that a larger 
number of independent K-K relations connect the real and the imag- 
inary part of the function under examination. We also compute the 
sum rules for the suitable moments of the real and imaginary parts, 
and observe that only considering higher powers of the susceptibility 
the correct vanishing sum rules are more precisely obeyed. All our re- 
sults are in fundamental agreement to recent theoretical findings. Sum 
rules providing explicit information about structural properties of the 
material seem to require wider spectral range. These constraints are 
expected to hold for any material and provide fundamental tests of 
self-consistency that any experimental or model generated data have 
to obey; similar tests of coherence can be performed for other nonlinear 
optical processes, e.g. pump-and-probe. Verification of K-K relations 
and sum rules constitute unavoidable benchmarks for any investiga- 
tion that addresses the nonlinear response of matter to radiation on a 
wide spectral range. 
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1 Introduction 



Kramers-Kronig (K-K) dispersion relations and sum rules ^|21 El E] have con- 
stituted for a long time fundamental tools in the investigation of light-matter 
interaction phenomena in condensed matter, gases, molecules, and liquids 
because they provide constraints able to check the self-consistency of exper- 
imental or model-generated data j3]. In particular these general properties 
allow to frame peculiar phenomena related to matter-matter or light-matter 
coupling that are very relevant at given frequencies, such as excitonic or po- 
laritonic effects in solidsjO], in the context of the interaction on the whole 
spectral range, showing that their dispersive and absorptive contributions 
are connected to all the other contributions in the rest of the spectrumjT]. 
Moreover by applying K-K relations it is possible to invert optical data, i.e. 
acquiring knowledge on dispersive phenomena by measurements over all the 
spectrum of absorptive phenomena or viceversa [71 . The conceptual foun- 
dations of such general properties, is the principle of causalityjEI] in the 
light-matter interaction, and the Titmarsch's theorem^J provides the con- 
nection between the mathematical properties of the functions describing the 
physics in the time and frequency domains. 

In spite of the ever increasing scientific and technological relevance of op- 
tical harmonic generation processes, and in general of nonlinear optical phe- 
nomena, relatively little attention has been paid to experimental investigation 
of the general properties of the corresponding nonlinear susceptibilities [T^ 
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Uni en en] , because the research has usually focused on achieving high reso- 
lution in both experimental data and theoretical calculations. Rapidly devel- 
oping technologies of tunable lasers allow to span larger and larger spectral 
ranges, so that the experimental investigations of frequency dependent non- 
linear optical properties of matter is becoming more and more possible. In 
the context of this new class of experiments, K-K relations and sum rules 
could provide information on whether or not a coherent, common picture of 
the nonlinear properties of the material under investigation is available. 

The first heuristic applications of K-K theory to nonlinear susceptibilities 
date back to the '60 [121 El (HI, while a more systematic study has begun 
essentially in the last decade. Some authors have preferentially introduced 
the K-K in the context of ab-initio or model calculations of materials prop- 
erties [inii2Diinii22i- 

Other authors have used a more general approach able to give theoretical 
foundations of dispersion theory for nonlinear optics [231 123 12S1 1211 12Z| • The 
instruments of complex analysis have allowed the definition of necessary and 
sufficient conditions for the applicability of K-K relations, which require the 
holomophicity of the susceptibility in the upper complex plane of the relevant 
frequency variable, given the form of the nonlinear susceptibility descriptive 
of the nonlinear phenomena under examination. It has also been shown that 
the asymptotic behavior of the nonlinear susceptibility, which can be ob- 
tained using the nonlinear Kubo-response function formalism [2H1, determines 
the number of independent pairs of K-K relations that hold simultaneously. 
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In the cases K-K do not hold, other techniques can be apphed for the inver- 
sion of optical data [3 123- Combining K-K relations and knowledge of the 
asymptotic behavior of the nonlinear susceptibility sum rules for nonlinear 
optics, obtained by other authors by following different strategies (SU], can be 
naturally derived along the same lines of the linear case. 

Recent theoretical advances [SH E21 have generalized the results obtained 
for the second [2Zj and third harmonics [221 by showing arbitrary order 
harmonic generation susceptibilities ^^'"■'(ria'; w, ■ ■ ■ ,u) are holomorphic in 
the upper complex u plane, and asymptotically decrease as there- 
fore K-K relations hold for all the moments uj'^°'x^'^\nijj] co, . . . ,uj) (from now 
on instead of x'^"^(na;; to, . . . ,uj) we use the simpler notation x^^K^^)) with 
< a < n. From these a total of 2n+2 sum rules for the moments of the 
real and imaginary part of the susceptibilities can be derived, only one (the 
(2n+l)*^ of the imaginary part) being different from zero and descriptive of 
the inner structure of the material under investigation. The previous results 
have been recently extended jHl] to the wider class of the moments of the 
/c*^ powers of harmonic generation susceptibility, which can be written as 
uj'^^\x^'^\nuj)Y . In this case K-K relations hold if < a < k{n -|- 1) — 1, 
and therefore a total of 2 A; (n -|- 1) — 2 sum rules can be derived, with the 
{2k{n+ 1) — 1)*'^ moment of imaginary part giving the only nonzero summa- 
tion over all the spectrum. The fundamental reason of this results is that if 
the harmonic generation susceptibility is holomophic in the upper complex 
ui plane, so are its powers. Higher powers of the susceptibility have faster 



5 



asymptotic decrease, so that the hmitations related to the finite range of real 

data are expected to be relaxed: this should be of particular benefit to the 

convergence of the sum rules, for which the consideration of the asymptotic 

behavior is more critical than for K-K relations jSHl 13 EH] • 

The consideration of anchor pointsjHZlEH] in the unknown part of x*-"-* {nu) [321 

is a very promising technique for improving data inversion. 

Until now, only in few investigations independent measurements of the real 

and imaginary parts of the harmonic generation susceptibilities have been 

performed for a relatively wide range El 1121 HH]; and consequently the 

verification of the coherence of measured data by check of the self-consistency 

of K-K relations is still of very limited use|42j. 

In this paper we present the first analysis of harmonic generation data 
where the full potential of the generalized nonlinear K-K relations and sum 
rules for harmonic generation susceptibilities is exploited. We consider two 
published sets of wide spectral range experimental data of third harmonic 
generation susceptibility on different polymers, the polysilane |12I (frequency 
range: 0.4 - 2.5 eV) and the polythiophene [lOj (frequency range: 0.5 - 2.0 
eV). 

In section 2 we apply K-K transformations to the real and imaginary part 
of the moments of the third harmonic generation susceptibility uj'^^x'^'iuj), 
with a suitable choices of a, and describe the quality of the data inversion 
obtained in this way. We also perform the same analysis for a'^"[x^(3a')]^, 
in order to prove for the first time on experimental data that K-K relations 
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hold also for the powers of the susceptibility, and show how many additional 
independent dispersion relations connect the real and the imaginary part. 
In section 3 we compute the sum rules for the experimental data, by calculat- 
ing the integrals of the suitable moments of the real and imaginary parts of 
the susceptibihty and of the /c*^ power of the susceptibilities, with 1 < A; < 5, 
and present a discussion on the issue of their convergence. 
In section 4 we present our conclusions. 

2 Efficacy of generalized K-K relations for 
data inversion 

In this study we base our calculations on two published sets of experimental 
data on third harmonic generation on polymers, where the real and imagi- 
nary part of the susceptibility were independently measured. The first data 
set refers to measurements taken on polisylane jl2] and spans a frequency 
range of O.4 - 2.5 eV. The second data set refers to measurements taken on 
polythiophene [10] and spans a frequency range of 0.5 - 2.0 eV. 
In this paper we control the self-consistency of the two data sets by observing 
the efficacy of the K-K relations in inverting the optical data for the func- 
tions u;2°[x^(3u;; w,a;,cJ)]^ with k=l,2. The theoryjSHEaiSSlEl! prescribes 
convergence of the dispersion relations: 




(1) 
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only ifO<Q;<A;(r7, + l) — 1, therefore in our case of n = 3 K-K relations 
should work only if0<a;<3if/c = l and < a < 7 if we consider 
the second power of the susceptibility. In our study we do not assume any 
asymptotic behavior outside the data range, but use only the experimental 
data, because extrapolation is somewhat arbitrary and in K-K analysis can 
be quite problematic jHUHl; we effectively apply truncated K-K relations and 
use a self-consistent procedure. 

In the paper by Kishida et al. j32] a check of validity of the K-K relations 
was already performed by comparing measured and retrieved x*^^^ (Scj). Apart 
from the fact that our analysis considers also another data set, the considera- 
tion of the moments of the susceptibility is not a mere add-on to the work by 
Kishida et a/. |121, but it represents a fundamental conceptual improvement: 
these additional independent relations are peculiar to the nonlinear phenom- 
ena, and provide independent double-checks of the experimental data that 
must be obeyed in addition to the conventional K-K relations. 

We present in figures C] and El the results of K-K inversion for respectively 
the real and imaginary part of the third harmonic generation susceptibility 
data on polysilane. We observe that in both cases the retrieved data obtained 
with the choices a = 0, 1 are almost indistinguishable from the experimental 
data, while for a = 2 and a = 3, the agreement is quite poor in the lower 



8 



part of the spectrum: the error induced by the presence of the cut-off at the 
high frequency range becomes more critical in the data inversion for larger a, 
since a slower asymptotic decrease is realized. We expect that inverting the 
data with the additional information given by anchor points located in the 
lower part of the data range these divergences can be cured. However, here 
our object is to deal with the worst case i.e. there is no a priori information 
about the phase of the complex nonlinear susceptibility at one or more fixed 
anchor points (fixed angular frequencies) jSHl- From the theory we expect 
that for a = 4 no convergence should occur: actually we observe that, while 
the main features around 1.1 eV are represented, there is no convergence 
at all for lower frequencies; the absence of a clear transition in retrieving 
performance between the a = 3 and the a = 4 case is due to the finiteness 
of the data range. 

In figures El and 13 we show the comparison between retrieved and experi- 
mental data of third harmonic generation susceptibility on polythiophene for 
the real and imaginary part respectively. The dependence of the accuracy 
of quality of data inversion is similar to the previous case: for a = 0, 1 the 
agreement is virtually perfect, while for a = 2, 3 we have progressively worse 
performance in the low frequency range; anyway the peaks in the imaginary 
part are still well reproduced, while the dispersive structures in the real part 
are present but shifted towards lower values. In this case the quality of re- 
trieved data for a = 4 is more distinct from what obtained with a = 3 
than in the previous data set. The inversion with a = 4 presents a notable 
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disagreement in the whole lower half of the data range for both real and 
imaginary part, in particular we see that in figure El the dispersive structure 
is absent, while the main peak in figure E] is essentially missed. 

Usually it is likely to expect that only the real or the imaginary part of 
the nonlinear susceptibility has been measured. Then normal procedure is 
to try data inversion using K-K in order to calculate the missing part. 

The results on figures 1—4 confirm that best convergence is obtained 
when using conventional K-K, therefore these should generally be used to 
obtain a first best guess for the inversion of optical data, and should be 
used as seed for any self-consistent retrieval procedure; nevertheless if there 
is good agreement with the inversions obtained with higher values of a, it 
is reasonable to conclude that the dispersion relations provide much more 
robust results. In this sense, the two data sets here analyzed are quite good 
in terms of self-consistency. 

We underline that if on one side considering a higher power of the suscep- 
tibility implicitly filters out noise and errors in the tails of the data, on the 
other side experimental errors in the relevant features of the spectrum, peaks 
for the imaginary part and dispersive structures for the real part, are greatly 
enhanced if higher powers of the susceptibility are considered; in the latter 
case consistency between K-K inversion of different moments is expected to 
be more problematic than in the k = 1 case. Therefore improved convergence 
for more moments will occur for the powers of the susceptibility k > 1 only 
if the data are basically good. 
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In figures and ini we show the results of K-K inversion for respectively the 
real and imaginary part of the second power of the third harmonic generation 
susceptibility on polysilane. Up to our knowledge this is the first analysis 
of this kind on experimental data. We observe that for a = 0,1,2 the 
agreement between experimental and retrieved data is almost perfect, while 
it gets progressively worse for increasing a. Nevertheless as long as a < 6 the 
main features are well reproduced for both the real and imaginary part and 
the retrieved data match well if the photon energy is > 1.0 eV. The theory 
predicts convergence for a = 7 and divergence for a = 8: in our analysis we 
have divergence also in the former case, and it is reasonable to attribute this 
to the cut-off in the high frequency range of the data, because a very high 
moment as the M*'^ requires a very well defined asymptotic behavior. 

In figures [7| and |H1 we repeat the same analysis for the second power of 
the susceptibility data taken on the polythiophene: also in this case the 
agreement is very good if a = 0, 1, 2, but the narrower frequency range does 
not allow the data inversion if the very high moments are considered. If we 
consider the real part -figure 0- for a = 3,4 K-K data inversion provides a 
good reproduction of the experimental data for photon energies < 0.7 eV; 
for a > 5 there is no convergence in the lower half of the spectral range. For 
the imaginary part -figure |H1- we can repeat the same observations, except 
that for a = 5 there is still good reproduction of the main features of the 
curve. 

We emphasize that if experimentally only one of the real or imaginary 
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part of the harmonic generation susceptibihty has been measured, there is 
no direct use of K-K relations relative to higher powers of the susceptibility, 
since the multiplication mixes the real and the imaginary parts. Therefore 
in this case the K-K relations for A; > 1 can be used as tests of robustness of 
the results obtained with the dispersion relations applied to the conventional 
susceptibility. 

3 Verification of Sum Rules 

Sum rules for optical functions which obey K-K relations can be generally 
obtained by combining the dispersion relations with the knowledge of their 
asymptotic behavior, which can be obtained with a detailed analysis of the 
physics of the system at microscopic level, by applying the superconvergence 
theorem [TTl El |7j to the dispersion relations. 

In the case of arbitrary order harmonic generation processes, we have that 
for large values of angular frequency ~ ^puj~'^^~^ [S1IS21, where tp 

is a material-dependent constant. Therefore from the K-K relations (0) and 
Q it is possible to derive the following sum rules [5^ : 



oo 



w^"3f?{ [x^''^ {nu)]''}duj = 0, < a < A;(n + 1) - 1 



(3) 




(4) 
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poo 

/ a;2'=("+i)-i!3{[x^")(^^)]''}da; = (5) 
Jo 2 

The verification of linear sum rules [HI IH El Bj from experimental data is 
usually hard to obtain because of the critical contributions given by the out- 
of-range asymptotic part of the real or imaginary part of the susceptibility 
under examination P^I36j : however in the case of linear optics information of 
the response of the material to very high frequency radiation can be obtained 
with synchrotron radiation [2] . In general a good accuracy in the verification 
of sum rules is more difficult to achieve than for K-K relations, therefore a 
positive outcome of this test provide a very strong argument to support the 
quality and the coherence of the experimental data. 

In the case of harmonic nonlinear processes the technical limitations for 
achieving information for a very wide frequency range are very severe, and 
the verification of the sum rules is critical, especially for those involving 
relatively large values of a which determine a slower asymptotic decrease. 
Nevertheless, if we consider increasing values of k, the integrands in the 
equations Q, (Hj), and (0) have a much faster asymptotic decrease, so that 
the missing high-frequency tails tend to become negligible. Therefore we 
expect that for a given a the convergence of the sum rules should be more 
accurate for higher values of k, if we assume that the main features of the 
spectrum are well reproduced by the experimental data, as explained in the 
previous section. 
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We first focus on tlie vanisfiing sum rules In order to liave a 

measure of how precisely the vanishing sum rules are obeyed for the two 
experimental data sets under examination, we introduce the dimensionless 
quantities Zsr and Zcg-. 



' '/:r:c.2«|5R{[x(^)(3u;)]n|da;'' 



^ '/::;; a;-+i|3{[x(3)(3..)]n|d^'- 
Low values of Z^[a^ k) and Z<^{a, k) imply that the negative and positive 
contributions to the corresponding sum rule cancel out quite precisely com- 
pared to their total absolute value. The two data sets of the polymers have 
quite different performances in the verification of these sum rules. 

We present in figures 9 and 10 the results obtained with the data taken 
of polysilane by considering 1 < A; < 5 for respectively the sum rules of the 
real and the imaginary part: we can draw very similar conclusions for both 
cases. We see that generally for a given a, we have a better convergence 
when a higher k is considered, with a remarkable increase in the accuracy of 
the sum rules for A; > 3. Consistently with the argument that the speed of 
the asymptotic behavior is critical in determining the accuracy of the sum 
rule, we also have generally a decrease in the quality of the convergence to 
zero when, for a given k, one considers higher moments, thus increasing the 
value of a. Particularly impressive is the increase of performance in the 
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convergence of the sum rules of x^^K'^^) for both the real (2a = 0,2,4,6) 
and the imaginary part {2a + 1 = 1, 3, 5) when we consider k = 4,5 instead 
of = 1: the values of Zs^f and decrease of more than three orders of 
magnitude in all cases considered. 

In figures 11 and 12 we present the corresponding results for the experi- 
mental data taken of polythiophene. Most of the sum rules computed with 
this data set show a very poor convergence to zero, because the correspond- 
ing and ZsR are above 10~^; nevertheless we can draw conclusions similar 
to the previous case in terms of change of the accuracy of the convergence for 
different values of k and a; consistently with the relevance of the asymptotic 
behavior, the precision increases for increasing k and for decreasing a. But 
in this case for a given a the improvement in the convergence of the sum 
rules obtained by considering a high value of k instead of = 1 is generally 
small, being in most cases the decrease of and Zsr below or around an 
order of magnitude. 

We observe that the bias between the performances of the two data sets 
in the verification of the vanishing sum rules is extremely large, not com- 
parable to the discrepancies found in the analysis of the K-K relations. We 
can guess that the worse performance of the data on polythiophene can be 
mainly attributed to their less complete representation of the relevant non- 
linear electronic transitions of the material; as previously stated, the data 
extension is critical in the verification of sum rules. This result is consistent 
with the previously presented slightly worse performance of this data set in 
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the K-K inversion of the second power of x*^^)(3(jj), where the relevance of the 
out-of-range data is also quite prominent. 

Finally, from equation El it is possible to obtain a simple relation between 
the non vanishing sum rules for [x!"^\nu})Y and the A;*'^ power of the non 
vanishing sum rule for the conventional susceptibility x("^(?7,a'): 

For both of the two data sets considered in our work these relations do not 
hold for 1 < /c < 5. The equation |H1 relates the slowest converging sum rules 
for each A;, therefore it is reasonable to explain the poor performances of the 
experimental data in reproducing this theoretical results with the finite range 
of data under examination. 

4 Conclusions 

In this study we have performed the first thorough analysis of generalized 
K-K relations and sum rules on experimental data of third harmonic gener- 
ation. We have tested the consistency between the theory and experimental 
data in the worst case, which however is usually most typical, namely data 
on limited spectral range without extrapolations beyond the measured range, 
and without any knowledge of anchor points [HIl EHl IHnj- We have consid- 
ered two data sets of comparable spectral range referring to independent 
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measurements of 3fJ{x''^^(3ci;)} and Q'{x^^''(3a')} performed on two polymers, 
the polysilane whose data span a frequency range of O.4 - 2.5 eV and 
the polythiophene |1U], whose data span a frequency range of 0.5 - 2.0 eV 

We have inverted the optical data of the susceptibility using applying 
truncated K-K relations with a self-consistent procedure and have verified 
for the first time that K-K relations hold consistently also for the moments 
of the susceptibility 6ij^"x*^^)(3a;) with < a < 3, as predicted by the general 
theory [^ESIESI; the two data sets show very similar performances. The 
precision of the data retrieved with K-K relations is good in the upper part of 
the spectrum for all moments, while for a = 2,3 there is a disagreement with 
the experimental data in the very low portion of the spectrum, mainly due 
to the absence of data descriptive of the asymptotic behavior of the optical 
functions; nevertheless the main features of the data are well reproduced by 
the dispersion relations. The agreement between the results of the dispersion 
relations for the various moments of the susceptibility under examination 
provide fundamental information on the robustness of the retrieved data and 
on the self-consistency of the experimental data. 

We have repeated the same analysis for u'^"[x^^\3uj)]'^ , in order to give the 
first experimental confirmation of recent theoretical findings that predict that 
K-K relations hold for uj'^"[x^"'\nuj)]'' with < a< k(n+l) — l, so in our case 
dispersions relations connect the real and imaginary part of c<j2"[;^{3)(3t^)]2 
if < a < 7. We have found that for the experimental third harmonic 
susceptibility measured on polysilane K-K relations hold with a high degree 
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of precision for all the functions ci;^"[x''^''(3a;)]^ with < a < 6; the agreement 
between retrieved and experimental data is excellent for < a < 2. In the 
case of polythiophene data, we can present similar conclusions, except that 
the agreement exists only for < a < 5. The disagreement between the 
theory and the experimental fact can be safely attributed to the truncation 
occurring in the high frequency range, which mostly affects the convergence 
of the dispersion relations of the very high moments. K-K relations for higher 
powers of the susceptibility cannot be directly applied if measures on only 
one part of the susceptibility are available, but provide additional tests that 
inform on the completeness and self-consistency of measured and retrieved 
data. 

We have investigated the validity of the sum rules that can be derived 
by combining the knowledge of the asymptotic behavior of x*-^''(3u;) with the 
previously described sets of K-K relations for the moments of x'^^HS^) and 
of its powers. The theory predicts that the integration in the semi-infinite 
positive cij-domain of each even moments of the real part of the k^'^ power 
of x'-'"-'(?7,a') up to the {2k{n -|- 1) — 2)*^^ gives 0, and the same holds for the 
odd moments of the imaginary part of the k^^ power of the susceptibility up 
to the {2k{n -|- 1) — 3)*'^; the only nonvanishing sum rule is the given by the 
integration of the {2k{n -|- 1) — 1)*'^ moment of the imaginary part. In our 
case we have considered all the powers of x^'^\'^'^) to the fifth. 

In order to evaluate the performance of the sum rules that the theory 
predicts to be vanishing, we have introduced an dimensionless measure of 
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how precisely a sum rule computed on the experimental data converges to 
zero. Generally, for both the sum rules of the real and of the imaginary part, 
for a given a, the convergence improves if we consider a higher value of k. 
Moreover, a decrease in the precision of the sum rules is realized if for a given 
k, higher moments are considered. The main reason for this behavior is in 
the concept that the faster the asymptotic decrease of the integrand of the 
sum rule under examination, the smaller the error due to the high-frequency 
cut-off. The two data sets for polymers differ greatly in the precision achieved 
in the verification of the sum rules, in many correspondent cases the polysi- 
lane data provide results that are better by orders of magnitude. The main 
reason for this discrepancy, much more relevant but coherent with the results 
obtained in the K-K study of [x*-^'*(3a;)]^, is the much stronger dependence of 
the sum rules precision on the position of high frequency range experimental 
cut-off relative to the saturation of the electronic transitions of the material: 
it is likely that the data on polythiophene are, apart from being narrower in 
absolute terms, less comprehensive relatively to the electronic properties of 
the material. 

On the contrary, the experimental data on both polymers do not verify 
the non vanishing sum rules, which the theory predicts to give information 
about the structure of the material, because they do not obey the newly 
established consistency relation (jH)) which should hold for each value of the 
power k considered for the susceptibility. The non vanishing sum rules involve 
the highest moments for which the integrations converge, therefore they are 
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most strongly affected by the frequency range finiteness; improvements on 
the data range are then necessary to expect to obtain verification of these 
sum rules. 

The constraints here analyzed are in principle universal, since they es- 
sentially derive from the principle of causality in the response of the matter 
to the external radiation, and so are expected to hold for any material; they 
provide fundamental tests of self-consistency that any experimental or model 
generated data have to obey; similar tests of coherence can be performed for 
other nonlinear optical processes, e.g. pump-and-probe^l ITBl 1^ I44j . Verifi- 
cation of K-K relations and sum rules constitute unavoidable benchmarks for 
any investigation that addresses the nonlinear response of matter to radiation 
on a wide spectral range. 
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Figure captions 

Fig. 1: Efficacy of K-K relations in retrieving 3?{x*^^^(3a')} on polysilane 
Fig. 2: Efficacy of K-K relations in retrieving $5{x('^^(3a')} on polysilane 
Fig. 3: Efficacy of K-K relations in retrieving ^{x^^\'ioj)} on polythiophene 
Fig. 4: Efficacy of K-K relations in retrieving 53{x*^^''(3ci;)} on polythiophene 
Fig. 5: Efficacy of K-K relations in retrieving 3?[{x'-'^^(3a;)}]^ on polysilane 
Fig. 6: Efficacy of K-K relations in retrieving $5[{x^^^(3u;)}]^ on polysilane 
Fig. 7: Efficacy of K-K relations in retrieving 3f?[{x''^''(3u;)}]^ on polythio- 
phene 

Fig. 8: Efficacy of K-K relations in retrieving 53 [{x*-^-* (3"^)}]^ on polythio- 
phene 

Fig. 9: Convergence to of the vanishing sum rules w2"g?[{x(3)(3u;)}]'= with 
1 < A; < 5; data on polysilane 

Fig. 10: Convergence to of the vanishing sum rules a;^""'"^9'[{x'-^'*(3'^)}]^ 
with 1 < A; < 5; data on polysilane 

Fig. 11: Convergence to of the vanishing sum rules a'^°'3?[{x*^^)(3a')}]^ 
with 1 < k < 5; data on polythiophene 

Fig. 12: Convergence to of the vanishing sum rules ui'^"~^^^[{x^^\3u)}]^ 
with 1 < k < 5] data on polythiophene 
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Figure 2: Lucarini, The Journal of Chemical Physics 
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Figure 5: Lucarini, The Journal of Chemical Physics 
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Figure 6: Lucarini, The Journal of Chemical Physics 
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Figure 7: Lucarini, The Journal of Chemical Physics 
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